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Abstract. - We show that the asymmetric crystal environment of the V site in the ladder 
compound NaV2 05 leads to a strong coupling of vanadium 3d electrons to phonons. This 
coupling causes fluctuations of the charge on the V ions, and favors a transition to a charge- 
ordered state at low temperatures. In the low temperature phase the charge fluctuations 
modulate the spin-spin superexchange interaction, resulting in a strong spin-phonon coupling. 



Introduction. - The spin-Peierls transition in one-dimensional spin systems is one of the best 
C^ ' known and investigated cooperative phenomena in low-dimensional physics. This transition 

comprises a dimerization of the lattice of a Heisenberg spin chain, leading to an alternation of 
the exchange parameter J . As a result, nearest- neighbor spins form a singlet state accompanied 
by the opening of a spin gap. Consequently, the static spin susceptibility drops drastically 
below the transition temperature (Tsp)- This effect, predicted by Pytte [|l|, was first observed 
in one-dimensional organic compounds l2l . The discovery was followed by theoretical work of 
Cross and Fisher g , who showed a possibility to calculate Tsp from the coupling of spins to 
the lattice. The discovery of the spin-Peierls transition in the inorganic spin-chain compound 
CuGeOs B renewed interest in this phenomenon. In this half-filled chain compound each Cu 
site holds one "hole" with spin s = 1/2. The spin-spin exchange interaction occurs via oxygen 
orbitals. Lattice vibrations cause a modulation of the spin-spin superexchange due to changes 
in the inter- ionic distances and bond angles. This leads to a spin-lattice coupling which plays 
an important role as a driving force of the spin-Peierls transition. 

The ladder system NaV2 05 was proposed to be a second inorganic spin-Peierls compound 
[p|. It contains VO5 pyramids which form V-0 ladders extended along the 6— axis m, see 
Fig. la. For half the ladders the apical oxygen 03 is above the V ions while for the other half 
they are below them. From the chemical formula one finds an average vanadium valence of 
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+4.5. This implies that a 3d-electron is shared between two V ions within one rung, making 
NaV205 a "quarter-filled" ladder compound M. The shared electron introduces an additional 
degree of fi-eedom: there is a possibility for the electron to move between the right (R) and 
left {L) V ions within one rung. This feature leads to physical differences between NaV2 05 
and CuGeOa, the most striking one being a possible charge ordering in NaV2 05. If an 
additional degree of freedom is involved in a phase transition, one should expect enhancement 
of the anomalies accompanying it. In fact, available experimental data on NaV205 reveal 
strong anomalies in heat conductivity [g|, specific heat |^, inelastic neutron scattering [[lO| , 
and inelastic light scattering |l^. These observations can be considered as (i) confirmation 
of changes in the electronic subsystem in addition to a spin gap opening, and (ii) indirect 
evidence for charge ordering. Direct evidence of charge separation within a rung was provided 
by a nuclear magnetic resonance experiment |12| which showed that two non-equivalent V 
ions exist below Tc- In addition. X-ray and elastic neutron scattering revealed a lattice 
superstructure with a modulation vector (1/2,1/2,1/4) in the low-temperature phase [[l3| . 
However, the nature of the charge ordering and the strength of the charge asymmetry are not 
clear yet. 
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Fig. 1. - a) Schematic plot of a part the crystal structure of the NaV2 05 compound. Grey circles 
denote V ions while open circles correspond to the O ions. The V-O-V rungs are shown in bold, b) 
Charge disproportion arising due to out-of-phase V ion displacements along the z-axis. Light-grey 
(dark-grey) circle corresponds to decreased (increased) electronic density on the site. 



Two different models for the charge ordered phase were introduced recently, notably an 
"in-line" |1J] and a "zigzag" structure [1^, |l6|. In the in-line structure all V ions within one 
leg have the same charge Qa but Qil ^ QiB., where i enumerates rungs along the ladder, while 
I = L,R denotes the sites within one rung. For the zigzag state one has QiR — Qi+i.L- For 
both the states the electron moves within one rung, i.e. Qil + QiR = 1. It was proposed in 
[ [l4| that the phase transition occurs in two stages. First the charge subsystem orders, which 
is followed by the opening of a spin gap. Mostovoy and Khomskii IQ argued that a zigzag 
ordering itself may lead to the opening of a spin gap due to an induced alternation of the 
exchange within the ladders. Nevertheless, coupling to the lattice is always important for the 
spin-Peierls transition either as a driving force or as a factor stabilizing the charge ordering 
[ p6| . For a better understanding of the role the lattice plays it is important to know the values 
and physical origin of the parameters involved in a model taking the spin, charge and lattice 
degrees of freedom into account. 

Here we investigate a mechanism of electron-phonon coupling in NaV205 and determine 
the scale of the parameters of the spin-phonon interaction. Due to the asymmetric structural 
environment of the V sites (see Fig. la) the vibration-induced shifts of electron site energies 
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are proportional to z-axis displacements of the ions. The small V-03 equilibrium distance 
dap ~ 1-61 A m leads to a strong contribution of the 03 ion Coulomb field, E^p, to the 3d 
electron site energy. A phonon mode involving a non-equivalent modulation of dap on the right 
and left side therefore causes a charge redistribution between the V ions on a rung. It will 
be shown that in the charge-ordered phase this effect leads to a strong spin-phonon coupling 
activated by the charge disproportion. The Coulomb field enhances the role of electronic 
correlations proposed to be responsible for the charge ordering [^ Q |l6| and leads to 
phonon-driven fluctuations of the exchange. Finally, two possible experimental consequences 
of the charge disproportion and fluctuations for optical spectroscopy will be analyzed. 

One-rung charge disproportion. - Consider a V-O-V rung, denoting the indirect hopping 
between the V ions as t± « 0.35 eV [0, |l^. Then there are two in-rung states separated in 
energy by 2t±: an odd and an even one, and both V ions have charge -1-4.5. Now consider 
an out-of-phase vibration of the two V ions along the z-axis (see Fig. lb). The corresponding 
zero-point amplitude is zq ~ y^h/Mfl, where M is the V ion mass, and fl is the phonon 
frequency. For fl = 400 cm~^, which is an estimate justified later by our experimental data, 
zo ~ 0.04 A. Under a displacement (±z) of the V ions the on-site energies will be changed 
by an amount ±Cz, inducing a charge transfer between the V ions (C is the deformation 
potential). By determining the eigenstates of this asymmetric two-level system the ground 
state energy Sq and amount of the charge transfer is easily calculated, yielding 

1 Cz 



where Q = for V"*"^ state. A point-charge Coulomb estimate of -Eap = 2e/d^p yields C = 
gE^p = 12 eV/A. Therefore Czo/t±, the parameter characterizing the electron-phonon cou- 
pling, is of order unity. Assuming that an adiabatic approximation is valid {Q <^ t±), and calcu- 
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lating the mean value with the vibrational ground state wave function |0) = (vrz, 

we obtain the charge fiuctuations in the charge-symmetrical (high-temperature) phase: 



(01 (g - 1/2)^ \^) = \- ije«^^^" [1 - *(e/^o)] , (2) 

where ^ = tj_/C = 0.03 A, and $(^/zo) is the error function. Since ^ « zq, we obtain 

{(j)\ {Q — 1/2)^ 10) w 0.18, i.e. the fiuctuations are of the order of the net average charge. 

Now we can estimate the role of the displacement-induced site asymmetry in the formation of 

the low-temperature phase. Suppose that inter-site Coulomb interactions, leading to intersite 

correlations, induce a static splitting 2Acor of the on-site levels within a rung. This splitting 

causes a charge disproportion ip = \Ql — Qr.\- This purely correlation- induced disproportion 

ipcoT is of the order of Acor/i_L- However, the apical oxygen field at the d— electron considerably 

enhances this charge redistribution. The correlation-induced charge asymmetry will cause a 

difference in the Coulomb forces acting on the V ions. This difference induces static opposite 

displacements of the ions, leading to a further increase of the charge disproportion. Eventually, 

this process will stop due to elastic energy of the lattice. The additional splitting 2Aiat of the 

on-site levels leads to a total site energy shift ±A ~ Aiat -|- Acor, the average ion charges are 

1 /? 
given by Ql,r — 1/2 = ±A/2 (i^ -I- A^) . The new equilibrium positions z of the V ions are 

determined by a self-consistency equation: 

C^ 1 + . ^3) 



2tj_k y^l + T^{l + xy 
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where t = Acor/i±, ^ut = Cz, x = Aiat/ Acor, and k the non-d— electron contribution to the 
total force constant k — MQ^ (w 30 eV/A ) relevant for the vibration. Eq. (3) is a consequence 
of the equilibrium condition fcz = Cip/2 and Eq. (1) for the charge transfer between the V 
ions. It has the following limiting cases for r ^ 1: 

A/{l-A), 1-^2 ^ ^2 

X = ■{ , (4) 

21/3^-2/3 _i^ 1 _ ^2 <<- ^2 

with A = C'^/2t±k. The first limit corresponds to a response of the lattice to a relatively 
small correlation-induced charge disproportion. The second one shows that if the lattice is 
close to an instability (A fn 1) even weak correlations induce a large charge disproportion. 
The d-electron contribution to the force constant can be calculated from Eq. (1) as fc^; = 

d^eo/2dz^ « ~C^/2t^ (w -200 cV/A^ for C = eSap), yielding A = ~kd/{k ~ kj) ~ 0.85. 
It follows from Eq. (4) that x can be larger than unity, in which case the main contribution 
to (p comes from the apical oxygen (03) field. This field enhances the correlation induced 
asymmetry of the V ions and thus favors charge ordering. 

For T > 1 Eq. (3) yields x ^ 1/t. The charge disproportion and ion displacement now 
saturate at (^ « 1 and z ~ C/k, respectively. The corresponding z is of the order of zq, in 
agreement with the experimental shifts reported in Ref. [Q. The corresponding site energy 
shift is of the order of zC ~ t± . Within a linear approximation the phonon- induced charge 
transfer in the ordered phase can be written as: 

6Q - CSz ^^ ,,, , (5) 

where dz is the out-of-phase displacement from the new equilibrium positions. As we will see 
below, these charge fluctuations lead to a spin-phonon coupling in the low-temperature phase. 

Phonon-induced fluctuations of the exchange. - It has been proposed that above T^ the 
spin system of NaV205 can be mapped onto an antiferromagnetic Heisenberg chain with 
H = J^SiSi+i, where s^ is the spin per rung |0, and J ~ 560 K [g. Using this idea 
we can calculate J within a simple model. The superexchange occurs mainly due to hopping 
via the oxygen orbitals within the legs. We denote the hopping integral between V and O 
orbitals along the ladder as tpd, and the one-electron energy difference for the in- leg transition 
0(2p)— >V(3d) transition as e « 6.5 eV Q. Allowing for inequivalent nearest-neighbor rungs 
i and i + 1, J can be expressed as 

J.25:^^, (6) 

where the ti^i and ii+i,; are hoppings of a hole from the V sites to the O ion between them. 
The values of tf ^ being proportional to the probability to find the hole at a corresponding V 
site can be readily expressed via the charges of the ions: tf ^ = ti^Qi^i. From Eq. (6) we obtain 



J — '^tpdiQi,LQi+i,L + Qi,RQi+i,B.)/£'^t±{^. The charge-ordering modified exchange J for the 



(^) Expressing of J as a function of the charges is a common feature of the approaches based on 
perturbation theory M . 
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in-line and zigzag phases is uniform and may be written as 

T f ti, zigzag phase, Qj^L = Qi+i.ij 

J=jr-^-\ ■ (7) 

^ [ 2A^ + 1]^, in-Une phase, Qi^l = Qi+i,L 

The charge separation in the zigzag phase decreases the exchange while the in-line ordering 
increases it [^ . These effects can be easily understood since in the zigzag phase the electrons 
are more separated in space while in the in-line phase they are more close to each other. The 
experimentally found J/ J « 0.8 pO| thus is indicative of a zigzag type of charge ordering in 
NaVaOs. 

The phonon causing charge fluctuations also modulates the hopping tij and, in turn, the 
exchange J. For phonon momentum q along the ladder the modulation of J is found to be 

^-L=4J{QL~Qn)f{q)f-, (8) 

oz oz 

with f(q) = isin(qd/2) and /(g) = cos{qd/2) for the zigzag and in-line phases, respectively. 
SQ is the charge transfer in Eq. (5) and d the in-ladder lattice constant. This contribution 
to the spin-phonon coupling is solely due to charge asymmetry and therefore a characteristic 
feature of the low-tcmpcrature phase. Following the approach of Pytte Q| we estimate the 
spin-phonon coupling constant as A ^ Zq{S J / 6 z)'^ {nj)~^ . Note that, within the present model, 
a complete charge ordering (ip = 1,A^ t±) suppresses spin-phonon coupling since it prohibits 
phonon-induced charge fluctuations (see Eq. (5)). From Eqs. (5)-(8) we obtain: 

/^2 2 T A 2 

A^^^ ^ (9) 

The maximal A is reached at A = t±/\/2, corresponding to ip — l/VS, J/ J = 2/3, and A ~ 1. 
Finally we address the influence of charge ordering on possible frustrating next-nearest- 
neighbor interactions (J'). This is interesting since frustration {J' /J) may be a driving force 
for the opening of a spin gap. If we suppose that both nearest and next-nearest couplings 



are determined by the charge order, we obtain for the zigzag phase: J = J(l — ip'^) |19 
and J' — J' [(1 — ip)'^ + (1 -|- tp)^] /2 = J'(l + ip'^). Zigzag ordering apparently leads to an 
enhancement of the role of frustration. 

Experimental consequences. - The effects of the phonon-induced charge fluctuations and 
charge disproportion directly influence the electronic and lattice properties of the crystal and 
should therefore manifest themselves in Raman and infrared spectroscopy. 

One experimental consequence of the charge disproportion considered above is that the 
out-of-phase vibration of the V ions should become Raman active below Tc- The normal 
vectors of even (in-phase) and odd vibrations of the V ions (see Fig. lb) can be written as 
Ug = (1, l)/-\/2 and Uu = (1,— 1)/\/2, respectively, li ip ^ 0, these vibrations are mixed and 
their eigenvectors can be written in the form: 

Ui = Ug cos 6* -I- Uu sin 6*, U2 = — Ug sinff -I- Uu cos6'. (10) 

The angle 9 is determined by the charge disproportion within a rung [9 ~ at T > Tc). There 
are two reasons for Raman activity of the second vibration in Eq. (10). The first one is the 
presence of the Raman active Ug component in the U2 mode. The second reason is that charge 
asymmetry leads to electronically different V ions (Y+^-^+vf^ ^j^d V'*''*-^"'''/^). Therefore, 
displacements of these ions give different contributions to the Raman polarizability, leading to 
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Raman activity of the Uu mode. These two effects, which can not be separated presently, lead 
to an intensity I2 of the new mode proportional to ^ ip^ . 

Experimental Raman data for T = 100 K and T — 5 K are presented in Fig. 2. At 
T = 100 K, only a single Lorentzian mode at 420 cm^^ is observed. Below the phase 
transition, at T = 5 K, this mode splits in two components centered at 394 cm~^ and 429 cm~^, 
respectively. The lower energy peak appearing for T < Tc is attributed to the U2 vibration in 
Eq. (10). The ratio of the integrated intensities of the "new" (394 cm"!) and "old" (429 cm"!) 
phonons is I2/I1 ~ 0.4. This is considerably larger than for the phonons arising solely due the 
superstructure formation, where the ratio is found to be smaller than 0.1 ||2l|| . In principle the 
ratio I2/I1 could be used to estimate the extend of the charge ordering {ip ^ \fhji\)- However, 
the reliability of this estimate would be questionable since the Raman polarizability of the Uu 
is not only determined by charge ordering but is also influenced by resonance enhancement 
effects I2I. 




Raman Shift (cm" ) 



Fig. 2. - Raman spectra of NaV205 at T = 100 K (thin line) and T = 5 K (bold line) demonstrating 
the appearance of a new intensive phonon at 394 cm~^ in the low-temperature phase. Polarizations 
of the incident and scattered light are parallel to the ladders. The dashed line shows Lorentzian fits 
to the r = 5 K spectrum. 



Another experimental consequence is related to charge fluctuations as observed in infrared 
(IR) experiments. IR absorption spectra at T = 300 K show a peak at ?iWinax ~1 eV |g^. This 
peak has been assigned to in-rung electronic transitions. The in-rung hopping was estimated 
there as t\_ « 0.3 eV, whereas the difference ?iw,nax— 2t_L ~ 0.4 eV was interpreted as a 
result of an in-rung charge separation already at T > T^- Here, we propose that this shift 
may be due to the electron-phonon coupling considered above, even without static charge 
separation. In the presence of the out-of-phase phonon mode the splitting of the one- rung levels 
is 2y/t\ + C^z^ > 2tj^. Taking into account that the probability of a z— axis displacement is 
determined by the vibrational wave function \(f>), we obtain a band in the absorption spectra 
which starts at 2t± and peaks at ?iWmax ~ 2t^ + {(j)\C'^z'^/t± \(j)). Taking the same set of 
parameters as before, we get hujniax— 2tj^ w C'^z^/t^ w 0.6 eV, in good agreement with the 
experiment. Based on this model we expect an increase of TiWinax in the low-temperature 
phase, as well as a narrowing of the peak due to the suppression of charge fluctuations by the 
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charge ordering. 



Conclusion. - We investigated electron- and spin-phonon coupling in NaV2 05 caused 
by the structural vanadium site asymmetry. The electron-phonon coupling leads to charge 
fluctuations, and considerably increases a correlation induced charge disproportion. The 
fluctuations induce a strong spin-phonon coupling in the charge-ordered phase. We stress, 
however, that this interaction is not the driving force of the transition. The charge ordering 
at low temperatures is manifested in Raman scattering by a new intense phonon peak. The 
phonon- induced charge fluctuations in the high-temperature phase may explain the absorption 
peak observed near 1 eV in IR experiments. 
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